Density matrix expansion for the MDI interaction 
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By assuming that the isospin- and momentum-dependent MDI interaction has a form similar to 
the Gogny-like effective two-body interaction with a Yukawa finite-range term and the momentum 
dependence only originates from the finite-range exchange interaction, we determine its parameters 
by comparing the predicted potential energy density functional in uniform nuclear matter with 
what has been usually given and used extensively in transport models for studying isospin effects in 
intermediate-energy heavy-ion collisions as well as in investigating the properties of hot asymmetric 
nuclear matter and neutron star matter. We then use the density matrix expansion to derive from 
the resulting finite-range exchange interaction an effective Skyrme-like zero-range interaction with 
density-dependent parameters. As an application, we study the transition density and pressure at 
the inner edge of neutron star crusts using the stability conditions derived from the linearized Vlasov 
equation for the neutron star matter. 
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I. INTRODUCTION 

One of the phenomcnological nuclcon-nucleon inter- 
actions that have been extensively used in transport 
models to study heavy-ion collisions is the isospin- and 
momentum-dependent MDI interaction [l| . By using this 
interaction in the isospin-depcndcnt Boltzmann-Uhling- 
Uhlcnbcck (IBUU) transport model, an extensive amount 
of works have been carried out to study various isospin 
sensitive obscrvables in intermediate-energy heavy-ion 
collisions (for recent reviews see Refs. @>S|)- From com- 
parisons of the isospin diffusion data from NSCL-MSU 
for 124 Sn+ 112 Sn reactions at E = 50 MeV/nucleon Q 
with results from the IBUU model, a relatively strin- 
gent constraint on the density dependence of the nuclear 
symmetry energy at subsaturation densities has been ob- 
tained [f| [|| ■ The resulting symmetry energy has further 
been used to impose constraints on both the parameters 
in the Skyrme effective interactions and the neutron skin 
thickness of heavy nuclei Q • The MDI interaction with 
the constrained isospin dependence has also been used to 
study the properties of hot asymmetric nuclear matter [8[ 
as well as those of neutron stars [SHlIl]. including the 
transition density which separates their liquid core from 
their inner crust [15j . New constraints on the masses and 
radii of neutron stars were then obtained from comparing 
the resulting crustal fraction of the moment of inertia of 
neutron stars with that of the Vela pulsar extracted from 
its glitches [15j . 

Although the momentum-dependent nuclear mean- 
field potential or the energy density functional of uni- 
form nuclear matter was used in the above mentioned 
applications of the MDI interaction, the explicit form 
of the MDI interaction has never been given in the lit- 
eratures. It was, however, implicitly mentioned |l|, [l6[ 
that this momentum-dependent nucleon mean-field po- 
tential was constructed according to a Gogny-like inter- 
action |17| consisting of a zero-range Skyrme-like inter- 



action [18| and a finite-range Yukawa interaction as in 
the MDYI interaction [lg|. Similar to the Gogny inter- 
action, the momentum dependence in the nucleon mean- 
field potential from the MDI interaction is then from the 
contribution of the finite-range Yukawa interaction to the 
exchange energy density of uniform nuclear matter. The 
effect of the finite-range part of the MDI interaction in 
nonuniform nuclear matter was in previous applications 
either neglected or inconsistently included by using the 
average density-gradient terms from the phenomenolog- 
ical Skyrme interactions. From comparing the potential 
energy density functional obtained from above assumed 
interaction with that used in previous applications of the 
MDI interaction, we are able to find unique relations 
among the parameters of this interaction and those used 
in the potential energy density functional for uniform nu- 
clear matter. 

As a first step to facility the consistent application 
of the finite-range MDI interaction, it is of interest to 
derive an effective zero-range interaction with density- 
dependent parameters usingthe density matrix expan- 
sion first introduced in Ref . [19| . As shown in Ref. [2fJ for 
other finite-range interactions, corresponding zero-range 
interactions derived from the density matrix expansion 
reproduce reasonably well in the self-consistent Hartree- 
Fock approach the binding energies and radii of nuclei 
obtained with the original interactions, particularly if 
the density matrix expansion is only used for the ex- 
change energy from the finite-range interaction [19|, [2(| ■ 
In the present paper, we carry out such a study as an at- 
tempt to include consistently the effect due to the finite- 
range part of the MDI interaction. As an application, 
we use the resulting interaction to study the transition 
density and pressure in neutron stars and compare the 
results with those from previous studies using the MDI 
interaction but with inconsistent density-gradient coeffi- 
cients [H,[2l[. 

This paper is organized as follows. In Scc.[IIJ we review 



the isospin- and momentum-dependent MDI interaction 
and determine its underlying nucleon-nucleon (NN) in- 
teraction by fitting the energy density functional used in 
previous applications of the MDI interaction. In Sec. lIIIl 
the density matrix expansion is then used to derive from 
the finite-range exchange interaction a zero-range effec- 
tive interaction with density-dependent parameters. The 
application of the resulting interaction to study the tran- 
sition density and pressure at the inner edge of neutron 
star crusts is given in Sec. IIV1 using the stability condi- 
tions that are derived from the linearized Vlasov equa- 
tion for the neutron star matter. Finally, the summary is 
given in Sec. [Vj For details on the determination of the 
underlying nucleon-nucleon interaction potential in the 
MDI interaction and the application of the density ma- 
trix expansion to the exchange energy of nuclear matter, 
they are given in Appendices [X] and [Bl respectively. 



the finite-range exchange term in the nucleon-nucleon 
interaction, the intrinsic momentum dependence in the 
nucleon-nucleon interaction and the effects of short-range 
nucleon-nucleon correlations |25[ can also contribute to 
its momentum dependence. For simplicity, we assume in 
the present study that all of the momentum dependence 
in the MDI interaction comes from the finite-range ex- 
change term. In this case, the explicit form for the MDI 
interaction can be obtained from the energy density given 
in Eq. ([1]) by assuming that the interaction potential be- 
tween two nucleons located at r\ and ?i has a form simi- 
lar to the Gogny interaction [l|, [l7| but with its Gaussian 
form in the finite-range term replaced by a Yukawa form, 
that is 



v(fi,r 2 ) = h 3 (l + x 3 P a )p a f^±^) 



II. THE MDI INTERACTION 

In applications of the momentum-dependent MDI in- 
teraction for studying isospin effects in intermediate- 
energy heavy-ion collisions as well as the properties of 
hot asymmetric nuclear matter and neutron star matter, 
one usually uses following potential energy density 
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for infinite nuclear matter of density p and isospin asym- 
metry S — (p n ~p p )/p, with p n and p p being, respectively, 
the neutron and proton densities. In the above, t(t') is 
the nucleon isospin; f T (r,p) is the nucleon phase-space 
distribution function at position r; and po = 0.16 fm~ 3 is 
the saturation density of normal nuclear matter. Values 
of the parameters A\ = (A L + A u )/2, A 2 = (At - A u )/2, 
B, er, A, Ci = C T . r and C u = C r - T can be found in 
Refs. [l|, |5j. For symmetric nuclear matter, this inter- 
action gives a binding energy of —16 MeV per nucleon 
and an incompressibility Kq of 212 MeV at saturation 
density. 

The parameter x in Eq. (JIJ is used to model the 
density dependence of the symmetry energy. Based 
on analyses of the isospin diffusion in intermediate- 
energy heavy-ion collisions and the neutron skin thick- 
ness of heavy nuclei |^-|6J], the slope parameter L = 
3po(dEsym(p)/dp) p=Po of the symmetry energy has been 
constrained to L = 86 ± 25 MeV, which corresponds to 
— 1 < x < 0. More recent analyses of experimental data 
on isospin diffusion, double n/p ratio, and neutron skin 
thickness favor, however, a softer symmetry energy of 
L = 58 ±18 MeV [Hrll]. 

The momentum dependence in the MDI interaction 
can come from a number of different origins. Besides 
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In the above equation, the first term is the density- 
dependent zero-range interaction which can be consid- 
ered as an effective three-body interaction, whereas the 
second term is the density-independent finite-range in- 
teraction. In terms of this nucleon-nucleon interaction, 
the total potential energy of a nuclear system can be cal- 
culated from 



where 
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is the quantum state of nucleon i with the spatial state 
i r , the spin state i a , and the isospin state i T , and P r , 
P a and P T are, respectively, the space, spin and isospin 
exchange operators. 

As shown in Appendix^ where the detailed derivation 
of Eq. (PJ from Eq. @ is given, the first and second 
terms in Eq. ([1]) come from the direct contribution (the 
first term of Eq. ([3])) of the finite-range term in the NN 
interaction. The third term in Eq. ([T]) is simply from the 
contribution of the zero-range term. Although different 
density-dependent zero-range terms can lead to different 
symmetry energies/potentials [26J, we use in the present 
study the one in the standard Skyrme interaction [27j . 
The momentum-dependent terms in Eq. ([1]) come from 
the exchange contribution (the second term of Eq. §3§) 
of the finite-range interaction. Comparing the resulting 
energy density functional with Eq. (|T|) leads to following 
unique relations between the 8 parameters in the NN 
interaction and those in the energy density functional of 
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It is seen that the symmetry energy parameter x is re- 
lated to the coefficient of the spin-exchange term in the 
density-dependent zero-range interaction, whereas the 
width A in the momentum dependence is related to the 
range of the finite-range interaction. 

In the original MDI interaction PJ, the symmetry en- 
ergy parameter x can only have the value of or 1 . To get 
a larger range of density dependence for the symmetry 
energy while fixing the value of E sym (p ) = 30.5 MeV, 
the parameters Ai and A u have been expressed as [5| 
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which results in the x dependence of A%. Therefore, 
with this constraint the x dependence also appears in 
the finite-range direct interaction. 



III. THE DENSITY MATRIX EXPANSION 

Although the contribution of a finite-range interac- 
tion to the direct energy density of nuclear matter can 
be treated exactly, it is numerically challenging to eval- 
uate its contribution to the exchange energy density. 
The latter can be, however, approximated by that from 
a Skyrme-like zero-range interaction using the density- 
matrix expansion of Ref. [l9j . As shown in Appendix 
iBl the resulting exchange energy density from the finite- 
range interaction at position f in a nuclear system can 
be expressed in terms of densities p n and p p as well as 
kinetic energy densities r„ and t p as 

H§ L (r) = A[p n (r),p p {r)] 

+ B[p n (r),p p (r)]T n (r) + B[p p (f), p n (r)]T p (r) 

+ C{p n (r),p p (r)]\Vp n (r)\ 2 

+ C[p p (f), Pn (r)}\Wp p (f)\ 2 

+ D[p n (r),p p (r)]V Pn (r) ■ Vp p (r), (14) 
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B(p n ,p p ) = -2[p n V L (p n )+p p V u (p p ,p n )], 
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In above equations, Vnm(Pu, Pp) is the potential energy 
density of infinite nuclear matter from the finite-range ex- 
change interaction (i.e. the momentum-dependent terms 
in Eq. (JTJ) and is given by 



with 
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with ji and j'3 being, respectively, the first- and third- 
order spherical bessel functions and k T = (37r 2 p r ) 1 / 3 is 
the Fermi momentum. For the function F(p n ,p p ), it is 
defined as 

F(Pn,P P ) = \v L { P n)Pn + \v U {p p ,Pn)p p . (26) 



In Fig. Q] wc show the dependence of A(p n ,p p ), 
B(p n ,pp), C(p n ,p p ) and D(p n ,p p ) on p n and p p at 
subsaturation densities. It is seen that A(p n ,p p ) and 
D(p n ,p p ) are symmetric in p n and p p , while B(p n ,p p ) 
and C(p n ,p p ) are not. Furthermore, the density depen- 
dence of B(p n ,pp), C(p n ,Pp) and D(p n , p p ) are strong at 
low densities but weak near the saturation density. As 
these functions are from the exchange contribution of the 
finite-range interaction, they are independent of x. 




FIG. 1: (Color online) Dependence of A(p n ,p p ), B(p n ,p p ), 
C(p n ,Pp) and D(p n ,p p ) on p n and p p at subsaturation den- 
sities. 



IV. APPLICATION: TRANSITION DENSITY 
AND PRESSURE IN NEUTRON STARS 

The explicit form of the finite-range NN interaction 
and the Skyrme-like zero-range interaction with density- 
dependent parameters derived in previous sections make 
it convenient to use the MDI interaction for a wider range 
of studies. In this Section, we discuss its application 
in studying the transition density and pressure at the 
inner edge of neutron star crusts based on the stability 
conditions that are derived from the linearized Vlasov 
equation for the neutron star matter. 



A. The single-particle Hamiltonian 

The single-particle Hamiltonian for a nucleon can be 
obtained from minimizing the total energy of a nuclear 
system with respect to its wave function. For a neutron, 
it is given by 



K = -v 
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with the neutron effective mass 
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2m* 2m 
and the potential 

U n (p n ,p p ) = UZ(p n , P p)+U°(p n ,p p )+U*(p n ,p p ), (29) 
where 

dH (p ni Pp) 
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is the contribution from the the zero-range interaction 
with 

Ho(Pn,Pp) = ^t 3 p a [(2+x 3 )p 2 -(2x 3 + l)(p 2 n +p 2 p )] (31) 

being the corresponding energy density, and 

U°(p n , Pp ) = (w+*-H-M) fdh'p n {7)^^ 
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are the direct and exchange contributions from the finite- 
range interaction, respectively. Expressing the kinetic 
energy density in the exchange potential in terms of the 
density via the extended Thomas-Fermi (ETF) approxi- 
mation m, HI 



= apV z + b 



Pq 



cV 2 p q , q = n,p, 



(34) 



where a = |(37r 2 ) 2/3 , b = 1/36 and c = 1/3, the neutron 
potential can then be written as 

U n = U% + U° + UZ, (35) 

where 

K{Pn,Pp) = K{Pn,P P )+ dA{ ^ Pp) 

OPn 

+ a dB(p Pp ) /3 + a dB{p Pn ) /3; 

OPn Op n P 

UZ(p n , Pp ) = Gr\Vpn) 2 + G%>\Vp p ) 2 

+ G7V Pn ■ Vpp + G n n n2 V 2 Pn + G^ 2 V 2 Pp , 

(37) 



with 
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Similarly, the proton singlc-particlc Hamiltonian is 

h p = h„ (n^p) + U^ D + U% E , (43) 

where the first term h„ (n <-> p) is due to the nuclear in- 
teraction and is obtained from Eq. (|27|) by interchanging 
the neutron and proton symbols; the second term U^ 
is the direct Coulomb potential given by 

°?°-^ *'?T' ^ (44) 

and the last term UY E is the exchange Coulomb poten- 
tial which in the first order of the density matrix expan- 
sion [191 is given by 



U™ = -e 2 - 



1/3 



(45) 



B. The Vlasov equation for the neutron star 
matter 

The Vlasov equation has been widely used in studying 
the collective density fluctuation and spinodal instability 
in nuclear matter (Sec Ref. [3GJ for a recent review). For 
a /3-stable and charge neutral ripe matter, the Vlasov 
equation can be written as 



df q (r,P,t) 
dt 



+ v q ■ V ? f q (r,p, t) - V?U q ■ V?f g (r,p, t)=0 

(46) 



in terms of the Wigner function for particle type q 
n,p,e 
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where 4> q i is the wave function of i-th particle of type q. 
In Eq. (|4"6")l . v q — p/{m* + p 2 ) 1 / 2 denotes the particle 
velocity and to* = m e . 



To study the density fluctuation due to a collective 
mode with frequency w and wavevector k in the nuclear 
matter, we follow the standard procedure [30j by writing 



f q (r,p,t) = f° q (p) + Sf q (r,p,t) 



with 



8f g (r,P,t) =df q (p)e 
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As the momentum dependence of the Wigner function is 
mainly from f°{p), we have \/pf q (r,p,t) « Vpf°(f) and 
can thus rewrite the Vlasov equation (Eq. (14T>10 as 

- iu)Sf q + v q ■ {ik)5f q 



df° q (x-*U q 



where e q is the single-particle energy. Using Vpe q 
and writing p q {r, t) = p® + Sp q (r, t) with 
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which has same time and spatial dependence as 
Sf q (f,p,t), Eq. ([50]) can be rewritten as 



df q I^SUq 
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de 1 \ ^f 6 Pl' ' * / U-k-Vg 



(52) 



By substituting the above equation into Eq. (jSTj) , we ob- 
tain in the low-temperature limit 
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with s q = uj/kvg and t>^ = p E /(to* + p^ ) 1 ' 2 being 
the Fermi velocity. Carrying out the angular integration 
leads to the usual Lindhard function 
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In the low-temperature limit, the momentum integration 
can be approximately evaluated as 
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for q = n,p, with e q « pf? /2m* and 



X, 



Me 
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(56) 
pf is the electron chemical 



for electrons, where p, e 
potential. 

For the factor ^ , j-^Spqi in Eq. ([55)1 , there are the lo- 
cal contribution 5U®, the direct contribution SU^ from 
the finite-range interaction, and the gradient contribu- 
tion 5Un from the finite-range exchange interaction us- 
ing the density matrix expansion. In the case of neutrons, 
they are given by 
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after neglecting higher-order terms in 6p q . In the long- 
wavelength case, the finite-range direct contribution can 
be rewritten as 



sui. 
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DO 
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with 
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For protons, there are additional direct and exchange 
Coulomb contributions given, respectively, by 
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For electrons, there are only direct and exchange 
Coulomb contributions to ^ , -g-^-5p q >. 

After linearizing the Vlasov equation, the collective 
density fluctuation Sp q then satisfies the following equa- 
tion 



C f (Spn,Sp p ,Sp e ) T =0, 



(67) 



with 



C' 



XnLni-g^ + 

X T ( dU ° 4- 




eu 



DO 



dU,_ 

dp n dp 



9(7° ac/ D0 
) X p L p(~dpf + ~a^) " 




' ^„£n(Gr D + G™ 2 ) XnLniGg + Gf) 
X p L p (G™ D + G™ 2 ) X P L P (G^ D + Gf 2 ) 
On 




-Xe-L, 



4-7re v 
P~ A P 

^ 4?re 2 , 



U P 

SU? E 



\ 



— ) ) 



(68) 



r 



We note that the three terms in the above equation are 
due to, respectively, the bulk, the density-gradient, and 
the Coulomb contribution. 



C. The transition density in neutron stars 



Non-trivial solutions of Eq. (|67|) are obtained when 
\C'\ = 0, which also determines the dispersion relation 
u(k) of the collective density fluctuation. The transi- 



tion density in a neutron star is the density at which the 
collective density fluctuation would grow exponentially, 
resulting in the instability of the neutron star matter, 
and this happens when the frequency u> becomes imagi- 
nary. To determine the condition for this to occur, we let 
s q = —iv q (fg > 0) and rewrite the Lindhard function as 



Lq = -2 



2v q &rct&n(\/v q ). Since the values of L q 



are 
-2, 



—2 < Lq < 0, the critical values L n = L p = L e 
corresponding to v q — 0, then determine the spinodal 
boundary of the system when they are substituted into 



\Cf\ = 0. Expanding the nucleon and electron densities 
at low temperatures according to 
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we obtain from Eqs. (|55l) and (151)1) the following relations: 



l/2X e = fy*e/0p e , 

i/2X 9 = </a Pr 
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The bulk contribution in Eq. (|68[) in the low-temperature 
limit can thus be rewritten as 

X n (dp n /dp n ) X n (dp n /dp p ) 

8 ( X p (dnp/d Pn ) X p (dp p /d Pp ) |(73) 

X e {dp e /d Pe ) 



Denoting the effective density-gradient coefficients as 
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(75) 
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(77) 



we then obtain same expressions for determining the 
spinodal boundar y of rip e matter in the so-called dynam- 
ical approach [15], I3ll - l33| , if the Coulomb exchange terms 
are neglected. In particular, in the long-wavelength limit, 
the spinodal boundary is determined by the vanishing 
point of 



Vdyn(fc) = 

where 



\C f \L t =- 



X n X q X e {dp n /dp n 

dp p {dp p /dp n ) 2 



V + (3k 2 



4ne 2 



V = 

(3 -- 
( 



6U™ 



dp p 



n 



2 
TF 



dp n /dp n 8p p 

pp ~r yJ-^np ~r J-SpnjL) ~r UnnL, , 

dpp/dp n 
dp n /dp n ' 

47re 2 
(dp e /dp e + 6U^ E /S Pe )- 



k 2 ' 
(78) 



(79) 
(80) 
(81) 

(82) 



To find the upper limit of the spinodal boundary for all 
possible wavevectors, we first find the minimum value of 
Vdyn given by 



V dyn =V + 2(4 7 re 2 /3) 1 / 2 - (3k 



TFi 



(83) 



which occurs when k 2 = {Aire 2 / (3) 1 / 2 — k\ F . The density 
that makes Eq. (|83|) vanish then determines the spinodal 
boundary in the neutron star matter or the transition 
density at the inner edge of neutron star crusts. 
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FIG. 2: (Color online) Density dependence of effective 
density-gradient coefficients from different values of x for sym- 
metric nuclear matter. The finite-range direct and exchange 
contributions are compared, and the value extracted by Oya- 
matsu I34J] is also shown. 



In Fig. [21 we show the density dependence of effective 
density-gradient coefficients from different values of x for 
symmetric nuclear matter. It is seen that the contri- 
bution from the finite-range direct interaction is much 
larger than that from the finite-range exchange inter- 
action. Also, D nn decreases while D np increases with 



increasing value of x. Furthermore, D r 



D pp and 



D n p = D pn for symmetric nuclear matter, whereas all 
four D 's have different values for asymmetric nuclear 
matter. The density-gradient coefficients extracted from 
the MDI interaction arc, however, much larger than the 
value of 132 ± 12 MeV fm 5 shown in the figure that was 
used by Oyamatsu [34[ to fit the nuclear radii and is also 
consistent with the average value from different Skymre 
interactions. We note that the large density-gradient co- 
efficients in the presence study is due to our assumption 
that all of the momentum dependence in the MDI inter- 
action comes from the finite-range exchange term. 

In Fig. [3l we compare the transition density and pres- 
sure in this study with those from previous calcula- 
tions 111, H3 that use the value of 132 MeV fm 5 from 
Ref. [34[ for the density-gradient coefficients and neglect 
the exchange Coulomb interaction. It is seen that the 
present results are smaller than previous ones for neu- 
tron star temperatures T = MeV and T = 1 MeV. We 
note that including the exchange Coulomb term for pro- 
tons slightly increases pt while the larger density-gradient 
coefficients make pt smaller. From the latest constraint 
L = 58 ± 18 MeV on the density slope of the nuclear 
symmetry energy, the transition density and pressure 
are constrained within 0.050 fm~ 3 < p t < 0.071 fm~ 3 
and 0.12 MeV/fm 3 < P t < 0.31 MeV/fm 3 for T = 
MeV, and 0.038 fm~ 3 < p t < 0.070 fm" 3 and 0.06 
MeV/fm 3 < P t < 0.30 MeV/fm 3 for T = 1 MeV in the 
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FIG. 3: (Color online) The transition density p t and pressure 
Pt at the inner edge of neutron star crust as functions of the 
slope parameter L for T — MeV and T = 1 MeV. Previous 
results [lg, [21| based on density-gradient terms from Ref. [34| 
are also shown for comparison. 



term. As an application, the resulting interaction is used 
to determine the transition density and pressure at the 
inner edge of neutron star crusts based on the linearized 
Vlasov equation for the neutron star matter. Fortunately, 
the instability of neutron star matter is dominated by 
the bulk properties of nuclear matter, the large density- 
gradient coefficients due to the finite-range interaction do 
not reduce much the transition density and pressure in 
neutron stars. 



Appendix A: From the NN interaction to the 
potential energy density 



The two particle state \ij > in Eq. © contains the 
spatial, spin and isospin parts as 

\ij >= \iriairjrjajr >= \h-jr > \i<jj<j > \%t3t > ■ (Al) 

For the inner product of spatial state, we have 



present study. Although the transition density is smaller 
at fixed L compared to our previous results [la ], which 
leads to a smaller crustal fraction of the moment of in- 
ertia for neutron stars and an even stricter constraint on 
the masses and radii of neutron stars, the upper limit val- 
ues of p t and P t are larger because of the smaller values of 
L. The final constraint on the neutron star mass-radius 
(M-R) relation is expected to be similar to that in the 
previous work [15j . 



V. SUMMARY 

Assuming that the momentum dependence in the MDI 
interaction is entirely due to the finite-range exchange 
term in the nucleon-nucleon interaction, we have identi- 
fied the NN interaction potential that underlies the MDI 
interaction, which has been extensively used in study- 
ing isospin effects in intermediate-energy heavy-ion col- 
lisions as well as the properties of hot nuclear matter 
and neutron star matter through the resulting poten- 
tial energy density in infinite nuclear matter. Using 
the density matrix expansion, we have obtained an ef- 
fective zero-range Skyrmc-likc interaction with density- 
dependent parameters from the finite-range exchange in- 
teraction. Compared to the density-gradient coefficients 
extracted by Oyamatsu [34[ and the average value from 
different Skyrme interactions, the values from the MDI 
interaction are much larger. Since the momentum de- 
pendence in the MDI interaction could also come from 
the intrinsic momentum dependence in the elementary 
nucleon-nucleon interaction or be induced by the ef- 
fects of short-range nucleon-nucleon correlations, we have 
thus overestimated the effect of the finite-range exchange 



i,3 



E / < irjr\rir 2 X rir 2 \i r j r > d 3 r 1 d 3 r 2 



hJ 



E / tiPiWiVMPJMiV'Pndf 



'2 



p{ri)p{r 2 )d i r l d i r 2 



(A2) 



with (j>i(r) =< r\i r > being the space wave function for 
particle i and the density p(f) = J2t 44 (r)^i(^)) an d 



hi 

= ^ / < i r jr\fif 2 >< rif 2 \j r i r > d 3 r 1 d 3 r 2 
hj 



p{ri,r 2 )p{r 2 ,ri)d rid r 2 



(A3) 



with the off-diagonal density p(fi, r 2 ) — J2i 4>i {^1)^(^2) , 
which reduces to the density when f± = f 2 . 

For the zero-range term vq = 7^3(1 + 
x 3 Pa)p a {^ 2 -)S(f 1 - f 2 ) in the MDI interaction, 
its direct contribution to the energy can be calculated 
using the fact that P a = 1/2 when considering the inner 
product of spin state in spin saturated matter, and the 



result is 

i,0 



= §|d 3 rirf 3 r 2 p(f=i)p(f 2 )p° 



a fn+r2 



x 1 



d 3 rHP{r) 



S(n - f 2 ) 



(A4) 



with 



< = |(i + fy« 



(A5) 



Its exchange contribution 

E ° = \ E < *i N(~ p » -P a Pr)\i3 > (A6) 



'■./ 



can be evaluated by using {l+xzP a ){— P a P T ) = —P<tPt — 
x 3 P T and replacing P T by 8^ j T when considering the 
inner product of isospin state, and the result is 



Eh 



*3 

12 



d 6 r 1 d i r 2 p(r ll r 2 )p(r 2 ,r 1 )p° 



a fn+r 2 



.1 \ l + S 2 

: 5{n-r 2 ) [ -+x 3 J ~^— 



= Jd 3 rH E (r) 



with 






ff^ ^ _J± f ±J_~_ 1 ^"+2 



(A7) 



(A8) 



and <5 is the isospin asymmetry. The total contribution of 
the zero-range interaction to the potential energy density 
is thus 



H - H a +H a - — p 



1 - 



1 + 2x 3 rfj 



(A9) 



Comparing it with Eq. ([T]) leads to Eqs. © — ([7]). 

For the contribution from the finite-range interaction 
in Eq. ([2]) to the energy, we use in the exchange term 
the relation (W + BP a - HP T - MP a P T )(-P a P T ) = 
M + HP a - BP T - WP a P T to obtain the following total 
contribution to the potential energy 

E r = - d 3 r 1 d 3 r 2 {[p n (r 1 )p n (f 2 ) + p p {ri)p p {r 2 )} 

x \w + --H- Y )+2p n (? l )p p (r 2 )\W+j 
+ \p n (r l i,r2)pn(r f 2,r f i) + Pp(n,r 2 )p p (f 2 ,fi)} 

+ 2 p„(fi, f 2 )p p (f 2 , n) \M + - J | ' (A10) 



with the first and second terms being the direct contribu- 
tion, and the third and fourth terms being the exchange 
contribution. 

To obtain the energy density functional from the finite- 
range interaction, we introduce the coordinate transfor- 
mation 



f=(f 1 +r 2 )/2, s = ri-r 2 . 



(All) 



For the direct contribution, we use the approximation of 
infinite nuclear matter with constant density and make 
use of the integral J d 3 se~ fis / s = 4ir/p 2 to obtain 

E E = J d 3 rH E {r) (A12) 

with 

(A13) 
For the exchange contribution, we express the density 
matrix in terms of the Fourier transform of the Wigncr 
function, that is 



r + ^?- S -) 



Substituting Eq. (|A14|) into Eq. (|A10[) and evaluating the 
Fourier transform of the Yukawa interaction according to 



d 3 se- l{ P-P">- s 



s e~» s 4tt 



iti + tf-py/n 



, (A15) 



the exchange contribution of the finite-range interaction 
to the potential energy is then 

E E = f d 3 rH E {r), (A16) 



with 



H E {r) 



2tt 

~2 



d 3 pd 3 p' 



p" j 1 + (p - p*) 2 / fi 2 
x (fn(r,p)fn(r,p)+ f P (f,p)f p (f,p)) 



M+X-B-W 

2 2 



M 



H 



2f n (r,p)f p (r,j?) 



(A17) 



Comparing Eqs. (|A13|) and (|A17|) with Eq. (JTJ) then al- 
lows one to obtain Eqs. (j5j) — (fl~2")l for the parameters in 
the finite-range Yukawa potential in the MDI interaction. 



Appendix B: Density matrix expansion of the 
finite-range exchange interaction 

Here we follow exactly the same procedure for the den- 
sity matrix expansion introduced in Ref. [19( . To obtain 



10 



the energy density functional in Eq. (fl~4)) from the ex- 
change contribution in Eq. (jA10|) . we first use the coor- 
dinate transformation given by Eq. ()Allj) to express the 
density matrix as 



with. p SL (k T s) = -]^ji{k T s), g(k T s) = 2 (k^ s )» 33(Ks) and 

T T (r) = YA \V(j) Ti (r)\ 2 being the kinetic energy density. 
By making the approximation 



Pr(ri,r 2 ) 
= pAr + -,r--) 

i 



(Bl) 



with Vi(2) acting on the first (second) term on the right, 
and (j) T i is the spatial wave function of particle i with 
isospin r. The angular average over the direction of s is 
then 



1 
sinh[f • (Vi - V 2 )] 



s s 

dfl s p T (f+-,f- -) 



|-(Vx-V 2 , 
Using the Bessel-function expansion 

+00 

xy 



^^(r)^(r). (B2) 



£^M = ig(4n + 3)j 2 „ +1 (x)g„(y 2 ) (B3) 



n=0 



with J2n+i the (2n + l)th-order spherical Bessel function 
and 



,(y 2 ) 
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i " 
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(4n + 2-20!y 2{n ~ 
l\(2n + 1 - l)\(2n + 1 - 2l)\ 



, (B4) 



1 
4tt 



dAl s p T {r+-,r- -)p T >{r - -,r + -) 



~ Pr(r)p SL (k T s)p T/ (r)p SL (k T/ s) 
+ PT{r)p SL (k T s)g(k T ,s)s 2 

]v 2 p Tl (r)-T T ,{r) + h 2 T ,p T ,{r) 
+ PT'(r)PsL( k T's)g{k T s)s 2 



-V 2 p T (r)-T T (r) + -k 2 p T (r) 



and using the relation 



F{p n ,p p )y 2 p ri 



dF(p„,p p ) 



(Vp, ; 



dF(p n ,p p ) 
P P 



we obtain 



E* = J d 3 rHi L (f), 



(B6) 



V Pn ■ Vp p . (B7) 



(B8) 



where Hg L (f) is the Skymrc-like zero-range energy den- 
sity from the finite-range exchange interaction defined in 

Eq. mD. 



one obtains by identifying x ~ fc T s and y ~ (Vi — 
V2)/2fc r with /c r being the Fermi momentum and keep- 
ing up to the third-order Bessel function the following 
result: 



1 

4~ 



s s 

d£l s p T {r+-,r- -) 



PsL{k T s)p T {f) + g(k T s)s 2 
]\7 2 p T {r)-T T {r) + ^k 2 p T {r) 



(B5) 
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